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Abstract 



We present a description of very dense vortex lattices in highly anisotropic 
multilayers, for high fields parallel to the layers. We show that a magnetic 
field inversion can occur away from the center of a vortex, provided the 
layers are sufficiently far apart. 



The structure of vortices in layered superconductors presents a fascinating 
array of peculiar characteristics, depending on the direction of the external mag- 
netic field. Some of these have been studied quite some time ago, while others 
more recently |l], Q]. The case where the external field is parallel to the layers 
is an especially interesting one. It was shown numerically |§, for example, that 
in layered superconductors with inequivalent layers immersed in a magnetic field 
parallel to the layers, the field around a vortex can be inverted. The inversion can 
be achieved only when the in-plane penetration depth of the weakly supercon- 
ducting layers is much greater than that of the strongly superconducting layers. 
This work was based on the assumption of an essentially uniform order parameter 
amplitude along the layers. 

The present paper aims to show that the above mentioned field inversion for 
vortices parallel to the layers can be a really generic effect in highly anisotropic 
multilayers, where the superconducting structure undulates along the z-axis, as 
long as the period of undulation is large enough, and as long as the applied mag- 
netic field is high enough. We show in particular that if the period of undulation 
d of the superconducting structure is much larger than a penetration depth A 
along the z-axis, then the field of a vortex is inverted away from the center of the 
vortex, for sufficiently high external fields. 

What happens essentially is that the field h of each vortex can go quite fast 
along the z-axis from a nonzero value h(0) at the center of the vortex to the zero 
value of the Meissner state, provided the penetration depth A is small enough. 
In fact, for sufficiently large values of d/X the field h could go to zero within a 
distance much shorter than d, where d is the distance between layers. Indeed, 
it will do so for dense vortex lattices with d^> A, where there are vortex centers 
halfway between the layers, because the superconducting layers expel the mag- 
netic field, limiting it to the interlayer regions. It turns out now that the value 
h(0) of the field at the center of the vortex, which lies halfway between layers, 
is approximately equal to the value H of the external field, for large H. Ideally, 
it would be zero on the superconducting layer due to screening, if the field can 



vary rapidly enough. This rapid variation is ensured if d is much greater than 
A, unlike the case of the usual superconductors or of high temperature super- 
conductors, where the field changes very slowly on the scale of d. Thus, if d/X 
is large enough, the field will start out at a high value halfway between layers 
and drop abruptly within just a few A's. It will have therefore a large negative 
initial slope, which makes it shoot down very abruptly. It will thus overshoot 
past zero towards negative values, before the existence of superconductivity on 
the layer forces it to come back up towards zero. The field therefore will decrease 
very rapidly from a large h(0) to a negative value, then go to zero within the su- 
perconducting layer, and then rise again in the next interlayer region, where the 
external field can penetrate easily. A field inversion can occur therefore, if d/X 
and H are sufficiently large, because the large negative slope of h makes it drop 
down to negative values just before becoming zero around the superconducting 
layer. 

In order to demonstrate the above mentioned field inversion, we shall use 
a very general model that takes into account the existence of a nonzero order 
parameter between the layers 0. We shall assume in particular that the Gibbs 
free energy of the multilayer in an external field H parallel to the x-axis is 
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Here z is the direction normal to the layers. 
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We shall measure x, y, z in terms of d, the period of the superconducting 
structure along the z-axis. We render a(z) dimensionless by taking out a dimen- 
sionful constant a, so that a(z)/a= a(z). We measure \1/ 2 in units of a//3, the 
vector potential A in units of hc/2ed, the magnetic fields in units of hc/2ed 2 , and 
the energies in units of d 3 a 2 //3. Thus the Gibbs free energy density takes the 
dimensionless form 
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H - ij£ - AM 2 + ^(V x A - Hf] , (2) 

where v — % /2Mad 2 , T 2 = M/m and A 2 = mc 2 (3 /Wire 2 a. The penetration 
depth A should not be confused with the experimentally measured penetration 
depth. Note furthermore that A depends on the choice of the dimensionful con- 
stant a, which is unspecified so far, while X 2 /i>d 2 does not. 

In this paper we shall examine very dense vortex lattices. In such lattices the 
vortex cells are practically rectangular. Vortices exist in every interlayer spacing, 
but in a squeezed triangular formation, due to vortex repulsion (see Fig. 1). The 
vortices in the vortex lattice are parallel to the x-axis, because the external field 
is assumed to be along that axis. The order parameter \& becomes maximum 
on each layer, the distance between successive layers being d, but vortices will 
appear only every two layers along the z-axis, due to their mutual repulsion. 
Thus each rectangular vortex cell will have a length of 2L along the y-axis, and 
2d along the z-axis. As the field H increases, the vortex cells get more narrow, in 
other words L decreases, but the structure along the z-axis remains unchanged. 
Indeed, due to their mutual repulsion the distance between two vortices along 
the z-axis cannot be reduced below 2d. We have thus a triangular lattice, with 
vortices situated in every interlayer spacing. The spatial periodicity of the vortex 
lattice results in the quantization of the magnetic flux contained in each unit cell. 
We anticipate therefore that L oc 1/H. 

We expect each vortex to be axisymmetric very close to its center. The vector 
potential and the magnetic field vary typically over distances of the order of A 
along the z-axis, but of order TA along the y-axis. We shall work consistently in 
the limit d ^> A, and r — »■ oo. In other words, our results will hold for highly 
anisotropic multilayers, with the layers sufficiently far apart. 

Very close to the vortex center, the field is axisymmetric. Thus, in the vicinity 
of the center, we expect the lines of constant field to be ellipses centered at 



the vortex center, i.e. h = h(Jz 2 + y 2 /T 2 ). Since the magnetic field decays 
to zero pretty fast, within just a few A's along the z-axis, the various vortices 



do not overlap in the z direction. Note also that ot(z) ~ Oi(Jz 2 + y 2 /T 2 ) for 
\z\ ^> L/Y > \y\/Y. Thus, if L <C Y, a(z) is axisymmetric practically everywhere. 
It is thus reasonable to assume that the axisymmetry of the magnetic field will 
be approximately true over most of the cell, as long as Y is very large. We can 
easily verify that the ansatz 

sin0 . , s . s 

Ay = —fr-Mp), (3) 

A z = cos <f>A(p), (4) 

with A x = 0, can lead to such an axisymmetric single vortex field. The anisotropic 
polar coordinates p, on the y-z plane are defined through the equations 



p = Jz 2 + y 2 /Y 2 , (5) 



sin0 = -, (6) 

P 

cos0 = — , (7) 

Tp 

assuming that the origin is the center of the vortex. 

We shall therefore assume that this ansatz for the vector potential holds over 
the whole cell. Since a(z) is really a(p) for very anisotropic multilayers (Y — ► oo), 
provided \y\ <^C Y, the order parameter amplitude and the magnetic field are 
essentially functions of p only, justifying thus the ansatz of Eqs. (3) and (4). 

The magnetic field of the vortex will be along the x-axis: 

with 

Q(p) = -l+pA(p). (9) 

The approximate axisymmetry of the magnetic field and of a(z) imply that the 
magnitude of the order parameter is approximately axisymmetric as well: 

M/(p,0)=^(p)e^. (10) 



The Gibbs free energy density takes the following form, in terms of ip and Q: 

The superconducting structure is represented by a(z). If this structure has a 
period d, then a(z) is periodic with a period d. We expect that the vortices are 
situated halfway between the layers. Therefore a(z) will have an extremum at 
the center of the vortex. Combined with the periodicity, this implies that a(z) is 
an even function of z. 

We repeat that since r — »■ oo, we can assume that a(z) is really a(p), as long 
as \y\ < L <C T. Thus ip, Q, A and h are essentially functions of p. 

For a very dense vortex lattice, i.e. a high external field, the unit cell will 
be a truncated ellipse, with j/Cf within it (see Fig. 2), and it will resemble a 



rectangle with sides 2L and 2, where L <^T. In that case \z\ < p <\/z 2 + L 2 /T 2 k> 
\z\, hence p ~ \z\. The flux in the unit cell will then be 2L ■ 2 J hdp= 2tc. The 
field varies very little along the y-axis when the y coordinate is much less than T, 
so when L has become quite small the field h will be varying along the z-direction 
only, remaining almost constant along the y-direction. The order parameter will 
also be a function of z only, remaining practically constant in the y-direction. We 
can say equivalently that ip is a function of p only. 

The Gibbs free energy will be the integral of g over the area of the truncated 

cell of Fig. 2. Thus the Gibbs free energy per unit x-length will equal 

r L/r r i l 

/ dp2nTpg+ dp Tpg[2<n - 4co S - 1 (— )]. (12) 
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The equations that minimize this functional for < p < L/T are: 
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Indeed, these equations minimize the first integral in Eq. (12). Clearly, since 
L<r, these equations really hold only at the vicinity of the vortex center. We 
can easily show then that Q(0) = —1 and ip(0) = 0. Furthermore, dh/dp is zero 
at p = and ip is linear in p near the origin. 

Let us now look at the case L/Y < p < 1. If L < T, then the second 
integral of Eq. (23) can be approximated by the integral J L , r 4:Lg dp, because 
cos~ 1 (L/Yp)^i (tt/2) — (L/Yp) on most of the interval [L/Y, 1]. Thus we have to 
find the field equations that minimize J L , r g dp. These equations are 



, 3 t vQ 2 ip _ dhp 
p 2 dp 2 



a{p)i> + ^ + ^- = v-^, (15) 



where p « \z\ within the cell, and 



d (Vh-YHs d 2 Qij 2 

This last equation can be rewritten in the form 

pYf p -Yh + YH = ^ 2 . (17) 

We stress again that Eqs. (15) and (17) hold only in the interval L/Y < p < 1. 

Let us now discuss the boundary conditions. We already said that Q(0) = — 1 
at the center of the vortex. Also ip(0) = 0. Since the vortex is situated halfway 
between the layers, and since ip varies little along the y-axis, due to the small 
L and the large Y, the order parameter must be zero on all the planes that are 
located halfway between neighboring superconducting layers. But the edges of 
the cells lie also on such planes. Hence ^(1) = 0. The field h has to fit with that 
of the next cell, so we must also have that dh/dp be zero at the edge of the cell, 
i.e. h(l) = 0. Combined with the boundary condition if)(l) = and with Eq. 
(17) this implies that h(l) =H. The boundary conditions for Eqs. (15) and (17) 
are thus ^(1) = 0, h(l) =H at the edge of the cell. 

Since the interval [0, L/Y] is quite small, we can consider the boundary con- 
ditions ip(0) = 0, Q(0) = —1 as appropriate for Eqs. (15) and (17) as well, even 



though those equations apply strictly to the interval [L/T, 1] only. 

We conclude then that Eqs. (15) and (17) describe fully the dense vortex 
lattice if the y-length of the cell is much smaller than T, subject to the boundary 
conditions ^(0) = ^(1) = 0, Q(0) = -1, h(l) =H. 

We shall explore the qualitative consequences of this description through 
a simple variational model, before conducting a more careful numerical study. 
Thus, we shall adopt the trial order parameter 

where ipo is determined by the details of a(z), and we shall solve for Q. 

Equation (16), along with the boundary condition Q(0) = —1, yields 

Q(p) = -1 + THp 2 /2 + lV/3 (19) 

for < p < 1/2, as well as 

h(p) = H + fip (20) 

for 1/2 < p < 1, k and \x being integration constants. 

The boundary condition h(l) =H, i.e. h(l) = 0, gives p, = 0. Thus we get 

h(p) = H + k P , if < p< 1/2 

= H, if 1/2 < p < 1. (21) 

Let us now integrate Eq. (16) from | — e to | + e, e being a positive number 
tending to zero. Using the ansatz of Eq. (18) we get the boundary condition at 
p = 1/2: 

Th{\+)-Th{\-)=2uQ{\), (22) 

with u = cPtfo/X 2 . We combine now Eqs. (19), (21) and (22) to get 



Hence 



m l-) = TH + *1 = rg ' 12 -"» + 24 ", (24) 

y 2 ! 2 12 + 2w y ' 



while h{\-\-) =H. Hence h changes discontinuously across the superconducting 
layer. A similar change was presented in [2]. 

We now note that h{\— ) will be negative if and only if uo > 12 and TH > 
24u/(uj — 12). Hence this naive model shows that the field h of the vortex will be 
inverted if the external field H is sufficiently high, and if d/X is sufficiently large. 

We shall now obtain further results with the aid of numerical work done with 
the choice 

a(z) = l-J2ae- b{z - n -^ )2 . (25) 

n 

We shall assume that there is a vortex at z=0, and another at z=2. The super- 
conducting layers are at z=±l/2, z=±3/2, z=±5/2, etc. We have solved Eqs. 
(15) and (17) numerically along the z-axis, in the interval [0,1], assuming solu- 
tions ip{p) and Q(p) with p ?» \z\ and <fi p» it/2. Since T is large, our solutions 
are valid away from the z-axis as well, in regions with ?/ < T. The boundary 
conditions are, at the center of the vortex, ip(0) = 0, Q(0) = — 1, as mentioned 
earlier, while ^(l) = 0, h(l) =H at the edge of the cell. 

We repeat here that the boundary conditions at the origin are actually the 
boundary conditions arising from Eqs. (13) and (14), since Eqs. (15) and (17) 
are valid in [L/T,l] only. If L is much smaller than T though, we may consider 
them as appropriate for Eqs. (15) and (17). 

We find that for sufficiently high d/X, and for sufficiently large H, the field is 
inverted along a substantial interval, but definitely before |. This is happening 
because the slope of h is very large when d/X and H are large, and hence h drops 
so low that it can become negative. We see also that the field is practically zero 
along the superconducting layer, giving thus an almost ideal Meissner effect there. 
After crossing the layer, the field h rises again, till it reaches the value H. In this 
final interval beyond p = 1/2 we seem to have the normal state. The behaviour 
described here can be seen in Figs. 3, 4, 5, 6, 7. 



The curves in Figs. 3 and 4 were obtained for the input parameters T = 100 



H=2, v = 0.0001, a = 2.5, b = 10, d/X = 40, and they correspond to L = 2.2, 
L/Y = 0.022. The solution holds in [L/T,l] strictly, so h has been extended 
towards the origin in Figs. 3a, 5 and 7, using the fact that it has zero slope there. 
The curves in Figs. 5 and 6 were obtained for the input parameters Y = 100, 
H=2, v = 0.0001, a = 2.5, b = 10, d/X = 200, and they correspond to L = 2.8, 
L/Y = 0.028. Finally the curve in Fig. 7 has the input parameters Y = 50, H=15, 
v = 0.0001, a = 2.5, b = 10, d/X = 40, and it corresponds to L/Y = 0.006. 

Note that the field inversion and the situation described by these Figures, 
where the vortex at z=0 and the normal state at z=l are connected through a 
superconducting region, were also discussed recently || in the context of ordinary 
superconductors in high fields. It should be noted furthermore that the above 
Figures were also obtained by imposing the boundary conditions V'(O) = ^(1) = 0, 
Q(0) = —1 and Q(l) = q, finding the solution and the corresponding Gibbs free 
energy, and then finding which Q(l) gave the minimum Gibbs free energy. We 
confirmed thus that the boundary condition h(l) =H arises naturally from the 
minimization of the Gibbs free energy. 

As H goes up, Q(l) and h(l) will increase. Thus the normal state at the cell 
boundaries gradually proceeds inwards, until Q(l) becomes sufficiently large to 
destroy the superconductivity on the layers. 

It should be mentioned that the above conclusions do not depend on the 
details of a(z). They depend on the fact that a(z) is even and periodic, with 
extrema at and 1, and a layer at |. Nor do they involve the y-distance between 
the vortices, since Y is large. 

We also note that the inversion results from the combination of two factors: 
the large negative slope of h at the origin, due to the large values of d/X and 
H, and the expulsion of the field from the layer. It should be therefore quite 
independent of the boundary conditions at z — 1. Indeed the equation for h 
has an attractive fixed point at h — 0. This is shown clearly in Figs. 3, 5, 7. 
We see that the field becomes zero the moment it starts crossing into the layer, 



maintaining this value up to the point when it starts exiting the layer. We have 
thus a perfect Meissner effect inside the layer. The inversion occurs around the 
point where the field starts entering the layer. It decreases linearly near the 
origin, with a large negative slope. This large slope makes it overshoot past zero, 
before turning upwards again, towards the zero value of the Meissner state. 

We have checked that the inversion is independent of the boundary conditions 
at z = 1 by solving Eqs. (15) and (17) from z = to z = 1/2, subject to the 
boundary conditions Q(0) = — 1, ip(0) = 0, and dh/dp = 0, dip/dp = at 
p = 1/2. We choose r = 50, H = 50, v = 0.0001, a = 2.5, b = 10, d/X = 40. The 
result, plotted in Fig. 8, shows again the inversion and the attractive fixed point 
at h — 0. 

We have shown then that in highly anisotropic multilayers, and for fields 
parallel to the layers, the field is inverted away from the center of a vortex, as 
long as d ^> A, i.e. as long as the vector potential is able to perform the necessary 
variations within the interlayer space, and as long as H is large. This result is 
insensitive to the boundary conditions at the edge. 

The field inversion discussed here could be seen experimentally in highly 
anisotropic multilayers, provided the layers are very far apart (d ^> A), at high 
fields. Construction of such multilayers would be especially interesting because 
it would be one of the few instances where the magnetic field would have ample 
space for varying. In all superconductors examined so far, even the high T c ones, 
the penetration depth is much longer than the interlayer distance. 

A field inversion has already been presented in [3]. However it was derived 
using a discretized model. The variation of \& along the z-axis, which is essential 
for the inversion to occur, is introduced in that work by having inequivalent layers, 
with various constant order parameter amplitudes. In our work the variation of ^ 
arises naturally from a(z), which creates the undulating structure. Furthermore, 
the result of [3] was derived near H cl , and it does not assume that d > A. Thus 



the inversion in [3] occurs mostly along the y-axis, unlike our inversion. 

A field inversion has been also presented in [6]. It applies however to non- 
layered anisotropic materials, and to tilted fields. Furthermore, the A there is 
large, unlike our A. 

Our field inversion occurs between the vortex and the normal state. It might 
occur though, for sufficiently large d/X, even earlier, when the elliptical cells just 
start touching. Thus, as the ellipse begins its truncation along the y-axis, it may 
be preferable for the field at the edges of the elliptic cell to become negative. Then 
the cell could be truncated along the y-axis without changing its flux! This field 
inversion would make possible a drastic reduction of the cell, without affecting 
the flux of each cell. Indeed, the flux of the cell would increase if we were to 
cut out some of the cell on the right and on the left, since the field would be 
negative in those pieces. We would need then to cut more on the left and on the 
right, pieces with positive h now, in order to bring the flux down to its quantized 
value. Thus the cell would be truncated considerably without affecting its flux. 
If such a field inversion at the edges of the cell were to occur, then we would have 
neighboring vortices with an inverted field in between. In such a case though, the 
lattice would be an array of parallel chains of vortices, and each chain would have 
a vortex in every interlayer spacing. The existence of such chains would require 
an attraction between the vortices along the z-direction, an attraction that would 
be made possible precisely because of the field inversion. 
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Figure Captions 

Figure 1: The vortex lattice at very high fields. Each cell has sides 2L and 2. 

Figure 2: A truncated elliptical cell. The distances are always in units of d. 

Figure 3: a) The field h as a function of p, for H=2, T = 100, a = 2.5, b = 10, 
d/X = 40, v = 0.0001. b) Detail of this graph. 

Figure 4: The order parameter amplitude if) as a function of p, for the parameters of 
Fig. 3. 

Figure 5: The field h as a function of p, for H=2, T = 100, a = 2.5, b = 10, d/A = 200, 
v = 0.0001. 

Figure 6: The order parameter amplitude if) as a function of p, for the parameters of 
Fig. 5. 

Figure 7: The field /iasa function of p, for H=15, r = 50, a = 2.5, b = 10, d/X = 40, 
i/ = 0.0001. 

Figure 8: The field h as a function of p, for H=50, T = 50, a = 2.5, 6 = 10, d/X = 
40, i/ = 0.0001, obtained with the boundary conditions dip/dp = and 
dh/dp = at p = 1/2. 



